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Preface 


The  U.S.  Army  Battlescale  Forecast  Model  (BFM)  is  a  mesoscale  short-range  forecasting  tool 
designed  to  run  on  a  single-processor  workstation  in  a  data-restricted  tactical  environment. 
Operated  at  resolutions  between  2  and  10  km,  the  model  is  hydrostatic  and  can  remain 
computationally  stable  at  relatively  large  time  steps  due  to  its  use  of  alternating-direction  implicit 
finite  differencing.  In  addition,  the  model  assimilates  data  using  Newtonian  relaxation  to 
incorporate  both  observations  and  the  time-tendencies  of  forecast  variables  from  a  previously 
run,  larger-scale  numerical  model.  The  BFM  is  currently  being  used  by  the  U.S.  Army  to 
produce  real-time,  short-range  mesoscale  forecasts  that  can  be  used  either  to  provide  input  for 
tactical  weather  decision  aids  or  to  assist  in  more  precise  firing  solutions  in  artillery.  Mean 
absolute  error  statistics  derived  from  a  small  sample  study  show  that  the  model  can  successfully 
assimilate  forecast  data  from  a  global  or  regional  numerical  model,  while  also  providing 
additional  forecast  resolution  within  the  boundary  layer.  A  few  case  studies  provide  examples  of 
how  the  BFM  can  simulate  low-level  diurnal  circulations  forced  by  terrain  inhomogeneities  and 
surface  radiation/heating  gradients  when  low-level  synoptic  forcing  is  weak.  By  maintaining  a 
timely  and  accurate  four-dimensional  gridded  database  of  meteorological  information,  a 
battlefield  commander  can  receive  critical  guidance  that  can  assist  in  determining  the  appropriate 
courses  of  action  while  in  a  rapidly  changing  battlefield  weather  environment. 
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Summary 


Introduction 

The  U.S.  Army  Battlescale  Forecast  Model  (BFM)  is  a  short-range  numerical  weather  prediction 
(NWP)  tool  designed  specifically  for  tactical  operations,  and  for  communication  and  data 
restricted  environments.  Featuring  a  simplified  graphical  user  interface  (GUI),  operators  can  use 
the  model  with  minimal  numerical  modeling  expertise.  This  model  allows  the  Staff  Weather 
Officer  (SWO)  to  view  meteorological  output  results  in  a  variety  of  graphical  forms. 

Furthermore,  operators  can  generate  on  demand  specifically  tailored  Tactical  Decision  Aids 
(TDAs)  or  weather  hazard  products  for  battlefield  execution  and  missing  planning.  The  model 
has  been  modified  and  simplified  for  operational  utility  in  tactical  environments,  aiming  for  a 
compromise  between  model  sophistication  and  timeliness.  This  report  will  demonstrate  that  the 
BFM  is  capable  of  producing  reliable  and  timely  mesoscale,  short-range  forecast  guidance  under 
most  typical  atmospheric  conditions. 

This  report  describes  the  BFM  from  many  perspectives,  including  physics,  numerics,  tactical 
application,  and  mesoscale  forecast  ability.  Additionally,  specific  case  studies  illustrate  the 
model’s  ability  to  simulate  local  mesoscale  circulations,  which  are  often  missed  by  the  coarser- 
resolution  models  run  at  the  various  operational  modeling  centers. 

The  BFM  is  a  mesoscale  NWP  software  system,  designed  to  be  globally  relocatable,  as  well  as 
robust  and  efficient  in  a  tactical  operation  setting.  The  model  is  able  to  initialize  coarser- 
resolution  NWP  model  output  supplied  from  U.S.  Air  Force  and  U.S.  Navy  operational  centers. 
The  system  is  also  capable  of  a  four-dimensional  data  assimilation  (4DDA),  based  on  the  method 
of  Newtonian  relaxation,  or  nudging.  Output  from  the  BFM  is  post-processed  by  an  external 
algorithm  called  the  Atmospheric  Sounding  Program  (ASP),  and  then  is  stored  as  four¬ 
dimensional  gridded  data  in  a  relational  database.  Once  in  the  database,  the  meteorological 
forecast  grids  can  be  rapidly  accessed  by  application  software  critical  to  tactical  operation 
decisionmaking,  such  as  weather  hazards  algorithms,  weather  visualization  tools,  and  TDAs. 
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1.  Introduction 


The  Battlescale  Forecast  Model  (BFM)  is  a  tailored  version  of  the  Higher  Order  Turbulence 
Model  for  Atmospheric  Circulation,  developed  at  Los  Alamos  National  Laboratory,  New 
Mexico.  The  BFM  was  selected  by  the  U.S.  Army  for  its  ability  to  simulate  boundary  layer  flows 
in  complex  terrain  (Yamada,  1981),  to  assimilate  observational  and  larger-scale  numerical  model 
data  (Yamada  and  Bunker,  1989),  and  to  run  efficiently  on  low-cost  workstations  and  personal 
computers  (Knapp,  1998).  Boundary  layer  meteorological  phenomena  that  can  be  resolved  by 
the  model  fall  between  the  meso-P  and  meso-y  classification  of  Orlanski  (1975). 

Hydrostatic  equilibrium  and  the  Boussinesq  approximation  are  assumed,  and  detailed 
parameterizations  for  both  surface  radiation  (Kondrat’yev,  1977)  and  boundary  layer  turbulence 
(Mellor  and  Yamada,  1982)  exist.  The  model  is  also  capable  of  treating  vegetative  canopy 
effects  (Yamada,  1982)  and  includes  a  grid-nesting  option  (Yamada  and  Bunker,  1988),  although 
they  are  not  used  in  the  current  BFM.  The  technique  of  Newtonian  relaxation  (Hoke  and  Anthes, 
1976;  Yamada  and  Bunker,  1989)  is  used  to  assimilate  both  observational  data  and  forecast 
parameter  tendencies  (from  a  previously  run  global  or  regional  numerical  model)  into  the 
prognostic  BFM  forecast. 

The  U.S.  Army  Integrated  Meteorological  System  (IMETS)  and  the  Computer  Assisted  Artillery 
Meteorology  (CAAM)  software  suite  are  U.S.  Army  systems  that  have  adopted  the  BFM  for 
real-time  tactical  operations.  Due  to  the  highly  transient  nature  of  the  atmosphere,  short-range 
mesoscale  forecasting  is  essential  to  mission  planning  for  global  U.S.  Army  operations  that 
require  detailed  knowledge  of  both  the  current  and  future  three-dimensional  meteorological  state. 
Many  of  these  operations  are  anchored  close  to  the  earth’s  surface  and  involve  areas  of  interest 
no  more  than  500  by  500  km  in  size.  These  operations  may  be  airborne  (ballistics,  munitions, 
rotorcraft,  fixed  wing),  ground-based  (troop  and  vehicle  mobility,  supply  and  transport,  gunnery, 
surveillance,  soldier  heat  stress),  or  strategically  placed  for  readiness  against  nuclear,  chemical, 
and  biological  releases.  Given  that  time  scales  of  mesoscale  phenomena  operate  on  the  order  of 
minutes  to  hours,  numerical  forecasts  that  provide  high-frequency,  gridded  output  are  critical  in 
supplying  meteorological  guidance  for  short-range  U.S.  Army  mission  planning  at  below-corps 
echelons. 

This  report  describes  the  BFM  and  discusses  its  applications  to  the  U.S.  Army,  examining  the 
model’s  strengths  and  weaknesses  under  different  large-scale  flow  patterns  and  regional 
topographies.  Included  is  a  summary  of  the  model  configuration  (terrain  ingest,  objective 
analysis,  and  initialization),  a  discussion  of  the  application  of  Newtonian  relaxation,  a  small- 
sample  statistical  evaluation,  and  a  few  brief  case  studies.  Significant  features  of  the  BFM, 
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particularly  the  assimilation  of  global  model  foreeast  tendencies  and  the  capability  to  simulate 
diurnal  mesoseale  boundary  layer  cireulations,  are  shown.  All  references  to  the  BFM  in  this 
report  relate  to  the  version  eurrently  operational  on  the  IMETS,  with  the  exeeption  of  a 
somewhat  different  configuration  used  by  the  CAAM,  described  in  section  2.5. 


2.  Model  Configuration  and  Initialization 


2,1  Model  Configuration 

Operating  on  workstations  or  personal  computers  that  have  been  selected  as  U.S.  Army  common 
hardware,  the  BFM  is  designed  to  generate  up  to  a  24-h  foreeast  depending  upon  the  availability 
of  larger-scale  numerical  model  guidanee.  Before  initiating  a  foreeast,  a  user  must  speeify  the 
model  domain  grid  center,  grid  dimensions,  horizontal  grid  spacing,  and  initialization  base  time 
from  a  graphieal  user  interface  (GUI)  window  (figure  1).  Onee  exeeuted,  the  BFM  module 
queries  the  IMETS  meteorologieal  database  for  available  upper-air  radiosonde  and  numerical 
model  data  for  objective  analysis,  within  an  expanded  region  surrounding  the  model  domain 
(valid  for  the  initialization  time).  This  region  is  a  161  by  161  grid,  spaeed  at  the  horizontal 
resolution  selected  by  the  user  (for  instance,  a  10-km  resolution  generates  a  1600  by  1600  km 
square  domain);  however,  surface  observations  are  colleeted  only  from  within  the  designated 
BFM  model  domain.  Once  the  input  meteorologieal  data  has  been  retrieved,  model  terrain  data  is 
generated  and  model  execution  initiated.  Upon  eompletion  and  post  proeessing,  the  output  data 
is  inserted  into  a  relational  gridded  meteorological  database  (GMDB),  after  whieh  the  user  can 
select  from  a  variety  of  visual  display  GUIs.  At  this  point,  BFM  output  can  be  displayed  or 
plotted  in  a  number  of  ways  (eontours,  cross  sections,  vertical  point  soundings,  tabular  format, 
SkewT-FogP)  on  any  available  IMETS  map  backgrounds. 
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BFM  Execution  Program  v5.0 


□a 


File  Options 


Help 


Grid  Center  Location 


Latitude  I43.850N 


MGRS  i34TBP92215854 


Longitude  I018.415E 


Map  Select 


Grid  Dimensions 
Grid  size  =  500.0  x  500.0  km 


Layout  51x51  i=i 


Spacing  10  km  a 


Initialization  Date/Time  07Apr  12Z  i=i 


NOGAPS  Data 

UA  Met  Data 
SFC  Met  Data 


i07Apr  12Z 
08Apr  12Z  Q 


View 


View 


Peaci  Met 

Run  BFM 

Figure  1.  Sample  GUI  for  the  BFM  execution  used  on  IMETS. 

The  IMETS  uses  a  horizontal  BFM  resolution  of  10  km,  with  16  terrain-following  vertieal  levels 
(32  in  the  CAAM  version)  and  a  rigid  top  set  7000  m  above  the  highest  terrain  elevation  grid 
point.  A  log-linear  vertical  stagger  is  applied  so  that  greater  vertical  resolution  is  retained  near 
the  surface.  The  BFM  terrain  elevation  data  is  generated  from  the  Defense  Mapping  Agency’s 
Fevel  1  (3  arc  s)  Digital  Terrain  Elevation  Data  by  taking  the  four  nearest  neighbor  elevation 
posts  relative  to  the  grid-volume  center  and  applying  bilinear  interpolation.  To  reduce  the 
possibility  of  spurious  gravity-wave  generation  in  the  presence  of  steep  terrain  gradients,  the 
model  applies  a  computational  smoother  to  the  prognostic  equations  at  each  time  step,  by  using 
the  four  surrounding  grid-point  values  and  a  predefined  smoothing  constant  (Yamada  and 
Bunker,  1989).  The  BFM  also  assimilates  external,  large-scale  numerical  forecast  data  that  acts 
to  damp  smaller-scale  wavelengths  that  may  erroneously  be  introduced  by  the  terrain. 

2,2  Objective  Analysis 

For  the  IMETS,  meteorological  input  data  used  by  the  BFM  currently  flows  one  way  from  an  Air 
Force  Weather  Agency  (AFWA)  Operation  Weather  Squadron  (OWS),  through  a  Tactical-Very 
Small  Aperture  Terminal  (T-VSAT),  and  down  to  a  fielded,  weather  effects  workstation  as 
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depicted  in  figure  2.  Raw  upper-air  radiosonde  and  surface  observations,  along  with  gridded 
numerical  model  output  produced  by  the  Naval  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS),  are  currently  received.  Details  of  NOGAPS  can  be  found  in  Hogan  and 
Rosmond  (1991). 

When  entering  the  BFM,  all  collected  radiosonde  data  is  subjected  to  gross  and  hypsometric 
checks  to  assure  quality  control.  Then  an  objective  analysis  is  executed  to  produce 
meteorological  fields  for  both  initialization  and  time-dependent  data  assimilation  (using 
NOGAPS  background  forecasts).  First,  the  radiosonde  data  is  interpolated  on  predefined 
Cartesian  height  surfaces,  after  which  spatial  interpolation  to  the  BFM  grid-point  locations  is 
accomplished  using  an  inverse-distance  weighting  function.  Then  the  data  is  reinterpolated 
vertically  from  the  Cartesian  surfaces  to  the  model,  terrain-following  coordinate  surfaces  defined 
by  eq  1 ; 


where 


z*  =  H 


Z-Zg 

H-Zg 


z*  =  transformed  vertical  coordinate 

z  =  Cartesian  vertical  coordinate 

zg  =  ground  elevation  above  mean  sea  level 


H 

H 

For  simplicity,  H 


=  material  surface  top  of  the  model 
=  corresponding  height  in  Cartesian  coordinates 
is  defined  in  eq  2  as 


H  —  H  4"  Zg  max 


where  Zg  max  is  the  maximum  value  of  Zg  on  the  grid. 


(1) 


(2) 


NOGAPS  forecast  data  is  analyzed  in  a  similar  fashion,  with  the  exception  that  a  two-pass, 
successive  corrections  scheme  (Barnes,  1964)  is  used  to  horizontally  interpolate  from  NOGAPS 
to  the  BFM  grid-point  locations.  When  both  radiosonde  and  NOGAPS  data  are  available,  a 
similar  successive  corrections  approach  is  used  to  combine  the  analyzed  first-guess  background 
field  (NOGAPS)  with  the  radiosonde  observations.  Weighting  in  the  compositing  scheme  has 
been  modified  slightly  from  that  described  previously  in  Henmi  and  Dumais  (1998)  in  order  to 
allow  radiosonde  data  to  have  greater  influence  at  grid  points  nearest  to  the  observation  location. 
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Figure  2.  Data  flow  from  AFWA  OWS  hub  to  IMETS,  with  Mbps  rates  for  the  continental  U.S. 

2.3  Initialization 

The  BFM  execution  cycle  begins  3  h  prior  to  initialization  base  time  in  order  to  dynamically 
“spin  up”  the  model  fields.  Nudging  is  applied  to  the  analyzed  fields  at  all  interior  BFM  grid 
points,  and  at  the  start  of  the  preforecast  cycle,  a  simple  logarithmic  wind  profile  is  assumed  that 
satisfies  the  condition  of  mass  consistency.  Surface  temperature,  wind,  and  dewpoint 
observations  are  also  assimilated  into  the  lowest  model  levels,  using  a  form  of  observational 
nudging  and  a  weighting  scheme  somewhat  different  from  that  used  for  the  three-dimensional 
nudging  of  the  analyses.  This  weighting  is  based  on  both  model  resolution  and  grid-point 
distance  from  the  surface  observation  site  (Henmi  and  Dumais,  1998). 

2.4  Post-Processing 

Upon  forecast  completion,  the  BFM  produces  a  number  of  model-state  parameters,  such  as 
the  u  and  v  horizontal  wind  vector  components,  potential  temperature,  and  water  vapor  mixing 
ratio.  While  these  prognostic  parameters  are  useful  for  U.S.  Army  operations,  it  is  necessary  to 
place  more  specific  weather  information  into  the  GMDB  for  tactical  software,  such  as  the 
Integrated  Weather  Effects  Decision  Aid  (IWEDA).  The  IWEDA  uses  this  output  to  provide 
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detailed  information  in  terms  of  why,  when,  and  how  weather  affeets  weapon  systems  (as  well  as 
their  subsystems  and  eomponents)  and  operations  (Sauter,  1996).  However,  the  IWEDA  requires 
data  not  direetly  supported  by  the  BFM,  sueh  as  turbulenee  and  surfaee  visibility.  In  order  to 
meet  this  data  requirement  a  post-proeessing  module,  the  Atmospherie  Sounding  Program 
(ASP),  was  developed  to  read  the  BFM  output  data  and  derive  additional  meteorologieal 
parameters 

for  the  GMDB.  Most  of  these  derived  parameters  are  used  in  U.S.  Army  aviation,  sueh 
parameters  as  turbulenee,  icing,  cloud  heights  and  amounts,  surface  visibility,  fog,  and 
thunderstorm  probabilities. 

Operational  forecasting  decisions  must  be  made  quickly,  and  it  helps  to  apply  artificial 
intelligence  (AI)  techniques.  Although  the  ASP  is  not  truly  an  AI  program,  the  module  applies 
elements  of  statistical  methods,  conventional  computer  programming,  and  basic  meteorological 
principles.  After  diagnostic  computations  are  complete,  the  operator  should  use  the  “if-then” 
rules  in  the  code  to  check  the  consistency  of  the  forecasts. 

Since  the  post-processed  products  must  be  applicable  at  any  point  in  the  world,  the  modules  are 
designed  to  cover  weather  events  at  any  location,  regardless  of  terrain,  elevation,  season,  time  of 
day,  or  climate  zone.  For  example,  Knapp  (1996)  collected  2790  surface  observations  from  80 
continental  U.S.  locations  from  July  1994  to  April  1995  and  at  all  hours  of  the  day.  Using  the 
surface  elevation,  temperature,  dewpoint,  relative  humidity,  windspeed,  ceiling  height,  and 
precipitation,  Knapp  implemented  screening  regression  techniques  using  stepwise  procedures  to 
determine  the  predictor  values  for  his  equations.  Once  the  best-correlated  predictor  was  found, 
other  predictors  were  included  to  achieve  optimal  statistical  results.  For  instance,  eq  3  is 
formulated  using  observations  with  derived  ceilings: 


VISIBILITY=8.06+  (0.0003*SFC  ELEVATION) 

-(0.0456*RH)+(0.0058*CEILING  HEIGHT) 


(3) 


where 

VISIBILITY 
SEC  ELEVATION 
RH 

CEILING  HEIGHT 


integer  value  category  (categories  in  miles) 
surface  elevation  in  meters  above  mean  sea  level 
relative  humidity  in  percent 
feet  above  ground  level 


Eigure  3  shows  a  plot  of  the  ASP-produced  surface  visibility  over  the  Midwestern  United  States, 
based  on  a  9-h  BEM  forecast  valid  at  2100  universal  time  coordinated  (UTC),  25  Eebruary  1999. 
A  more  detailed  discussion  of  the  ASP  is  included  in  Passner  (1999). 
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System 


Figure  3.  ASP  post-processed  surface  visibility  field  using  the  9-h  BFM  forecast  valid  at  2100  UTC, 
25  February  1999. 


2.5  CAAM  BFM 

CAAM  is  a  software  system  designed  to  enhance  the  accuracy  of  meteorological  (Met)  input  to 
artillery  aiming  calculations.  All  effective  artillery  firings  must  include  aiming  adjustments  that 
are  based  upon  Met  information  to  compensate  for  variations  in  atmospheric  wind,  temperature, 
and  density  along  the  artillery  shell  trajectory.  A  recent  study  of  the  component  errors  for  the 
Extended-Range  155-mm  artillery  round  concluded  that  for  a  trajectory  length  of  28  km, 

42  percent  of  the  range  error  (down-range  along  the  initial  firing  azimuth)  and  56  percent  of  the 
deflection  error  (cross-range  perpendicular  to  the  initial  firing  azimuth)  was  attributable  to 
inaccuracies  in  the  Met  data  (Reichelderfer  and  Barker,  1993).  Most  of  the  remaining  firing 
errors  are  due  to  muzzle  velocity  and  occasion-to-occasion  inaccuracies.  Therefore,  reducing  the 
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Met  errors  in  the  data  used  for  the  tiring  ealeulations  ean  signifieantly  improve  the  aiming 
accuracy  of  the  artillery  rounds. 

The  CAAM  BFM  was  used  to  generate  forecast  computer  meteorological  messages  (FCMMs) 
that  were  then  compared  against  those  derived  from  radiosonde  computer  meteorological 
messages  (RCMMs).  Because  the  CAAM  BFM  forecasts  data  rather  than  measures  it,  the  model 
predicts  Met  conditions  at  desired  points  in  space  and  time  and  has  the  potential  of  more 
accurately  representing  the  atmosphere,  thus  allowing  for  more  precise  artillery  aiming.  Since 
radiosonde  data  can  contain  inherent  inaccuracies  (such  as  balloon  drift)  for  artillery  aiming 
purposes,  an  alternate  approach  was  employed  for  comparison  that  makes  no  initial  assumption 
about  the  accuracy  of  the  data.  Met  analysis  was  performed  on  data  collected  during  Sense  and 
Destroy  Armor  (SAD ARM)  artillery  live  firings  that  occurred  at  Yuma  Proving  Ground,  AZ 
(YPG).  Actual  submunition  impact  data  were  compared  against  predicted  impacts  derived  from  a 
trajectory  simulation  program.  Two  types  of  Met  data  were  input  to  the  trajectory  simulator 
(RCMMs  and  FCMMs),  in  order  to  test  which  type  of  data  most  accurately  represented  the 
“real”  atmosphere. 

A  series  of  SAD  ARM  artillery  firings  were  simulated  using  General  Trajectory  Model  3 
(GTRAJ3)  simulations  from  the  Armament  Research  Development  and  Engineering  Centers 
Firing  Tables  Branch.  The  simulations  were  run  first  using  a  RCMM,  and  then  using  a  FCMM. 
The  specific  input  conditions  for  a  particular  firing  (SAD ARM  ballistic  and  aerodynamic 
parameters,  gun  azimuth/elevation  angle,  measured  muzzle)  were  held  identical  from  each 
GTRAJ3  run  and  were  compared  to  the  actual  impact  submunition  coordinates  obtained  at  the 
SAD  ARM  Reliability  Determination  and  Assurance  Program  (RDAP)  and  the  Limited  User  Test 
(LUT)  conducted  at  YPG  during  January,  April,  and  May  2000.  The  SAD  ARM  rounds 
experienced  and  responded  to  the  “true”  atmospheric  conditions.  Consequently,  the  GTRAJ3 
simulation  that  resulted  in  the  closest  impact  to  the  actual  impact  location(s)  was  the  one  using 
Met  data  most  representative  of  the  “true”  atmosphere.  The  RCMMs  were  approximately  3-h  old 
during  the  RDAP  (i.e.,  the  radiosonde  used  for  the  Met  message  to  aim  the  gun  was  taken  about 
3  h  before  the  firing).  The  FCMMs  were  3-h  forecasts,  valid  at  the  time  of  the  firing.  During  the 
LUT,  the  RCMMs  were  only  about  1.5-h  old  and  the  LCMMs  were  5-h  forecasts. 

Ligure  4  shows  the  actual  and  simulated  impact  points  for  one  of  the  guns  participating  in  the 
final  LUT  mission.  The  SAD  ARM  rounds  came  in  from  the  northeast  (from  the  upper  right 
toward  the  lower  left  in  the  plot,  as  indicated  by  the  diagonal  arrow).  Each  of  the  rounds  (labeled 
“Actuals”)  impacted  within  the  confines  of  the  target  array.  The  “+”  symbol  indicates  where  the 
GTRAJ3  simulation  hit,  using  a  1.5  h-old  RCMM  and  gun  I’s  aiming  data.  When  GTRAJ3  was 
rerun  using  identical  aiming  data  but  substituting  the  LCMM  (a  5-h  forecast  RCMM  that  was 
valid  at  the  firing  time),  the  impact  is  depicted  by  the  diamond  symbol.  As  shown,  the  GTRAJ3 
simulation  using  Met  data  from  the  CAAM  BLM  impacted  closer  to  the  actual  live  rounds, 
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indicating  that  the  CAAM  BFM,  even  when  forecasting  5  h  into  the  future,  more  accurately 
represented  the  true  atmospheric  conditions  along  the  trajectories  at  the  time  of  the  firing  than 
did  a  radiosonde  only  1.5-h  old. 
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Figure  4.  Actual  and  simulated  SAD  ARM  impact  points. 

Many  of  the  other  RDAP  and  LUX  firings  produced  similar  results.  Each  SAD  ARM  round 
contains  2  submunition  canisters,  thus  the  12  rounds  fired  during  the  RDAP  resulted  in  24 
ground  impacts  against  which  the  simulated  impacts  were  compared.  On  each  of  four  LUX 
mission  days,  24  rounds  were  fired,  resulting  192  total  impacts.  In  some  cases,  it  was  difficult  to 
visually  assess  which  simulation  was  closest  to  the  live  impacts.  In  order  to  quantify  the  results, 
the  straight-line  distance,  or  radial  miss  distance,  was  determined  from  every  submunition 
impact  coordinate  to  the  RCMM-GXRAJ3  and  FCMM-GXRAJ3  simulated  impact  points.  Mean 
radial  miss  distance  (MRMD)  values  were  then  calculated  for  the  two  RDAP  mission  days  and 
the  four  LUX  mission  days.  Xable  1  summarizes  the  MRMD  results  (values  are  in  meters): 

Table  1.  MRMD  summary  for  the  RDAP  and  LUT  SAD  ARM  missions. 


RCMM 

FCMM 

RDAP-1 

256 

222 

RDAP-2 

123 

48 

LUT-1 

235 

188 

LUT-2 

261 

306 

LUT-3 

178 

178 

LUT-4 

223 

141 
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As  per  table  1,  the  FCMM -based  MRMD  values  were  smaller  for  4  of  the  6  mission  days, 
indicating  that  the  FCMMs  more  accurately  represented  the  atmosphere  than  did  the  RCMMs  for 
a  majority  of  the  missions.  The  results  were  almost  identical  for  LUT-3  (when  rounded  to  the 
nearest  meter).  During  LUT-2  (18  Apr  2000),  a  disturbance  in  the  upper  atmosphere  caused 
strong  winds  (exceeding  50  ms'^  near  the  apogee  of  the  SAD  ARM  rounds)  that  changed  direction 
significantly  shortly  before  the  firing.  In  this  case,  the  RCMM  (1.5-h  old)  was  a  more  accurate 
representation  of  the  firing  conditions  than  the  5-h  FCMM  forecast. 


3.  Newtonian  Relaxation  and  Data  Assimilation 


Newtonian  relaxation  (nudging)  is  accomplished  by  introducing  an  additional  non-physical  term 
to  the  prognostic  equation  for  a  variable  a,  such  as  in  eq  4,  that  acts  to  “steer”  the  mesoscale 
model  solution  towards  an  observation(s)  or  large-scale  model  forecast  for  a  (Yamada  and 
Bunker,  1988;  Hoke  and  Anthes,  1976;  Harms  et  ah,  1992).  The  model’s  physical  forcing  terms 
are  included  in  F,  and  ao  is  the  observation  or  large-scale  numerical  model  value.  Although,  in 
the  case  of  winds,  ao  in  the  BFM  is  actually  a  target  wind  discussed  by  Yamada  and  Bunker 
(1989).  The  coefficient  Ga  is  used  to  determine  the  strength  of  the  nudging  (Pielke,  1984),  and 
this  value  in  mesoscale  models  is  somewhat  arbitrary  and  is  a  function  of  the  synoptic  situation, 
meteorological  parameters,  and  the  magnitude  of  other  model  forcing  terms.  Stauffer  and 
Seaman  (1990)  select  Ga  so  that  the  time  scale  of  the  slowest  physical  adjustment  process  in  the 
model  and  the  nudging  term  are  similar,  and  MESO  (1993)  shows  that  a  value  of  Ga  =  3  x  10'"^ 
implies  an  exponential  decay  towards  observation  with  an  e-folding  time  of  0.93  h  when  no 
model  physical  forcing  terms  are  present. 

da !  di  =  F(a,x,t)  +  Ga  *  (ao  -  a)  (4) 

For  the  BFM,  the  magnitude  of  Ga  is  dynamic  and  is  set  strong  enough  to  effectively  suppress 
undesirable  phase  and  amplitude  errors  and  correctly  assimilate  NOGAPS  fields  in  space  and 
time.  Conversely,  the  nudging  is  not  so  strong  as  to  significantly  inhibit  realistic  boundary  layer 
mesoscale  circulations  from  developing  within  the  model.  Values  of  Ga  in  the  BFM  range  from 
1  X  10'"^  to  3  X  10'"^,  although  when  nudging  surface  observations,  larger  values  are  used  due  to 
the  stronger  contributions  from  the  forcing  terms  in  the  prognostic  equations  (care  must  be  taken 
when  nudging  surface  temperature,  otherwise  the  model’s  simulated  boundary  layer 
development  may  be  stunted).  Values  of  Ga  in  the  BFM  are  determined  separately  for  the 
prognostic  variables  representing  wind,  perturbation  virtual  potential  temperature,  and  total 
water  vapor  mixing  ratio.  These  values  are  periodically  adjusted  within  the  BFM  (for  each 
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parameter)  during  the  model  foreeast  cyele  and  are  functions  of  the  height  above  ground,  time  of 
day,  BFM  boundary  layer  lapse  rate  forecast,  and  degree  of  synoptic  forcing  forecast  by 
NOGAPS.  For  example,  an  estimation  of  the  degree  of  synoptic  forcing  is  made  using  such 
identifiers  as  the  low-level  windspeed  and  the  surface  to  5000-magl  temperature  advection. 

In  conditions  of  weak  forcing  where  the  terrain  can  be  expected  to  have  a  strong  local  influence, 
the  magnitude  of  Ga  is  generally  reduced.  Local  flow  fields  forced  by  surface  temperature 
gradients  over  sloped  surfaces  (drainage/upslope  wind)  or  along  land/water  boundaries  (lake/sea 
breeze)  can  be  well  simulated  at  2  to  10  km  resolution  when  synoptic  flow  is  weak  (Doran  and 
Horst,  1982;  McNider,  1982;  Mahrer  and  Pielke,  1976).  Many  BFM  successes  have  been 
documented  under  such  conditions  (Yamada,  1981,  1983;  Yamada  and  Bunker,  1989;  Pielke  and 
Pearce,  1994). 

Three-dimensional,  hourly  gridded  “nudging”  files  are  created  by  linear  interpolation  between 
successive  6-h  NOGAPS  forecasts.  Although  the  spatial  resolution  ratio  between  the  1° 
NOGAPS  and  the  BFM  is  poor,  it  is  an  improvement  over  the  previously  used  2.5°-resolution 
NOGAPS  grids.  A  study  of  the  present  BFM  version  versus  a  non-hydrostatic  mesoscale  model 
used  for  operational  purposes  (Henmi,  2000)  indicates  that  comparable  and  often  improved 
surface  forecasts  can  be  obtained  by  the  “nested”  BFM. 


4.  Numerics  and  Physics 


A  typical  24-h  BFM  forecast  executes  in  45  to  60  min  on  a  single  processor  350  MHz 
workstation  (51x51x16  grid  points).  The  model  is  computationally  fast  due  to  the  use  of  an 
alternating  direction  implicit  (ADI)  solution  technique  (Peaceman  and  Rachford,  1955)  in  three- 
dimensions  (allowing  for  larger  time  steps),  in  addition  to  a  lack  of  explicit  moist  microphysics 
and  a  cumulus  parameterization.  The  ADI  method  is  second-order  accurate  for  both  space  and 
time  derivatives  and  is  unconditionally  stable  (Kao  and  Yamada,  1989).  Nevertheless,  special 
treatment  of  updating  values  at  intermediate  time  steps  is  required  to  apply  this  technique  in 
three  dimensions  (Richtmeyer  and  Morton,  1967).  By  using  the  ADI  method,  the  model  is  able 
to  reduce  the  speed  of  fast  external  gravity  waves,  significantly  increasing  the  maximum  time 
step  allowable  under  the  Courant-Friedrich-Lewy  (CFL)  numerical  stability  condition.  A  planar 
Universe  Transverse  Mercator  grid  is  used  along  with  a  staggered  C  grid  (Arakawa  and  Lamb, 
1977)  for  the  prognostic  calculations. 

The  pressure  gradient  in  the  planetary  boundary  layer  (PBL)  is  computed  by  first  assuming  that 
the  geostrophic  winds  at  the  PBL  top  (assumed  1500  magi)  and  at  levels  above  are  equivalent  to 
the  NOGAPS  forecasted  wind  values.  Below  the  PBL  top,  the  geostrophic  winds  are  computed 
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using  a  simplified  thermal  wind  relationship  similar  to  the  approximations  used  by  Yamada 
(1981).  An  advantage  of  this  method  is  that  realistie  boundary  layer  pressure  gradients  ean  be 
ealculated  with  less  likelihood  of  large  eomputational  errors  due  to  steep  model  terrain. 
However,  in  the  free  atmosphere  this  method  retards  the  possibility  of  the  BFM  generating  its 
own  pressure  gradient.  This  is  eurrently  not  so  problematic  since  relatively  strong  nudging  is 
applied  and  the  current  BFM  lacks  both  explicit  moist  physics  and  a  cumulus  parameterization. 

The  vertical  wind  component  is  determined  diagnostically  by  assuming  a  condition  of  mass 
continuity  and  no  vertical  motion  at  both  the  top  of  the  model  domain  and  the  earth’s  surface. 
Finally,  the  model  has  complete  parameterizations  for  surface  radiation  (Kondrat’yev,  1977), 
turbulence  (Mellor  and  Yamada,  1982),  stratus  cloud  radiation  (Hanson  and  Derr,  1987),  and 
stable  precipitation  (Sundqvist  et  ah,  1989).  A  more  complete  description  of  the  model  physics 
and  numerics  can  be  found  in  Henmi  and  Dumais  (1998)  or  Yamada  and  Bunker  (1989). 


5.  Forecast  Examples  and  Statistics 


5.1  Radiationally  Driven  Flows  Under  Weak  Synoptic  Forcing 

The  following  subsection  discussed  an  example  of  a  BFM  forecast  made  during  the  summer  in  a 
weakly  forced  synoptic  environment.  This  case  was  selected  to  illustrate  how  the  BFM  can 
simulate  boundary  layer  flows  that  are  primarily  forced  by  strong  surface  radiational  heating 
differences  (such  as  sea  breezes,  nocturnal  drainage  flows,  etc.).  Also  included  is  a  brief  mention 
of  a  second  forecast  example  produced  in  the  vicinity  of  the  Great  Salt  Lake,  UT. 

5.1.1  Case  Study  of  Florida  Sea  Breeze 

For  this  12-h  simulation  of  a  sea  breeze  circulation,  a  100  by  100  km  grid  of  2-km,  horizontal 
resolution  was  selected;  centered  near  Cape  Canaveral,  FL  (28.47°N,  80.53°W),  and  initialized 
at  1200  UTC  on  23  July  1999.  Since  the  synoptic  forcing  on  this  day  was  weak  and  the  large- 
scale  environmental  conditions  invariant  during  the  study  period,  archived  NOGAPS  model 
output  was  not  used  for  analysis  nudging — model  prognostic  equations  were  instead  nudged 
throughout  the  simulation  to  the  initial  fields  using  a  uniform  nudging  coefficient  value  of  3  x 
10'"^.  As  seen  by  the  1200  UTC  station  XMR  (Cape  Canaveral  radiosonde)  sounding  in  figure  5, 
the  ambient  low-level  pressure  gradient  was  quite  weak  with  southwesterly  flow  below  950  hPa 
and  easterly  flow  above  950  hPa,  with  environmental  windspeed  values  less  than  5  m/s  below 
550  hPa. 
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Figure  5.  Cape  Canaveral  sounding  data  valid  at  1200  UTC,  23  June  1999. 


Surface  values  for  albedo  and  soil  moisture  were  input  manually  based  on  values  listed  in  Pielke 
(1984)  for  various  types  of  surfaces.  Land/water  values  for  eaeh  grid  volume  were  also  determined 
manually  from  maps  of  the  area  with  a  diserete  categorization  of  either  100  percent  land  or  water. 
Water  temperatures  used  were  based  on  values  given  by  Rao  et  al.  (1999),  with  24  °C  used  for  the 
Atlantic  Ocean  and  26  °C  for  the  inland  water  bodies  (such  as  the  Indian  and  Banana  Rivers). 

Although  not  shown,  the  initial  surfaee  windspeeds  were  below  2.5  m/s  and  generally  out  of  a 
southwesterly  direction.  By  7  h  into  the  model  simulation  (figure  6),  a  well-developed  low-level 
sea  breeze  eireulation  had  developed  along  the  immediate  eoastline,  with  a  localized  secondary 
circulation  over  the  Indian  River  also  evident,  consistent  with  summer  season  observations  in 
this  area  (Atkins  and  Wakimoto,  1997)  and  with  eireulations  simulated  numerically  and 
discussed  in  Rao  et  al.  (1999).  Surface  (10  magi)  windspeeds  and  direetions  simulated  by  the 
BFM  appear  to  eompare  favorably  to  those  produced  by  the  ETA  Data  Analysis  System  (ED AS) 
analysis  (figure  7)  and  from  observations  provided  by  the  Kennedy  Space  Center  Mesonet 
(figure  8)  on  the  afternoon  on  23  July  1999. 
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Figure  6.  Forecast  (7  h)  of  surface  winds  (tn/s)  near  Cape  Canaveral,  valid  at  1900  UTC,  23  July  1999.  BFM 
resolution  is  2  km.  Contours  of  terrain  elevation  every  1  m. 
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Figure  7.  EDAS  surface  winds  valid  at  0000  UTC,  24  July  1999. 
Maximum  vector  is  9  m/s. 
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Figure  8.  Surface  wind  observations  (m/s)  from  the  Kennedy  Space  Center  Mesonet  valid 
at  1900  UTC,  23  July  1999.  Contours  of  terrain  at  every  6  m. 


5,1,2  Case  Study  of  Great  Salt  Lake  Breeze 

A  similar  12-h  model  run  was  produced  over  the  Great  Salt  Lake  at  2-km  resolution,  initialized 
at  1200  UTC  on  27  May  1999,  another  weak  synoptic  flow  day  (figure  9).  The  BFM  forecast  (7 
h)  shown  in  figure  10  captures  the  generation  of  a  well-defined  lake  breeze  front,  as  well  as  local 
upslope  circulations  due  to  the  sloping  orography.  These  flows  generally  compare  favorably  to 
corresponding  MesoWest  observations  shown  in  figure  1 1 . 
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Figure  9.  Salt  Lake  City  radiosonde  sounding  data  valid  at  1200  UTC,  27  May  1999. 
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Figure  10.  Forecast  (7  h)  of  surface  winds  (m/s)  near  Salt  Lake  City,  valid  at  1900  UTC,  27  May  1999.  BFM 
resolution  is  2  km.  Contours  of  temin  elevation  every  200  m. 
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Figure  11.  Surface  wind  observations  (m/s)  obtained  from  Utah  MesoWest  valid  at  1900  UTC, 
27  May  1999.  Contours  of  terrain  elevation  every  250  m. 


5,2  Statistical  Results  over  the  United  States:  Continental  U,S,  Study  from 
27  to  31  March  1998 

The  purpose  of  the  study  was  to  examine  the  BFM  performance  for  a  variety  of  geographical 
locations  during  an  active  synoptic  period.  The  following  13  pressure  levels  of  NOGAPS  data 
were  used  to  create  analysis  fields  for  nudging:  1000,  975,  925,  900,  850,  700,  500,  400,  300, 
250,  200,  150,  and  100  hPa.  A  total  of  7  BFM  forecasts  were  made  for  9  locations  within  the 
continental  U.S.  from  27  to  31  March  1998,  with  a  fixed  domain  size  of  500  by  500  km  and  at  a 
horizontal  resolution  of  10  km.  Each  forecast  was  24  h  in  duration,  and  a  statistical  analysis  of 
the  data  was  generated  every  6  h  for  surface  and  12  h  for  upper-air  parameters. 
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Plotted  ED  AS  analyses  of  the  500  hPa  height  and  wind  fields  are  shown  at  various  times  during 
the  study  period  in  figures  12  to  14,  which  illustrate  that  the  synoptic  pattern  during  the  period 
was  active.  Due  to  a  lack  of  surface  and  upper-air  observations  available  for  ground  truth,  the 
statistical  sample  includes  the  same  sites  used  in  the  BFM  initialization  (i.e.,  not  an  independent 
sample).  In  general,  about  3  to  5  radiosondes  and  8  to  15  surface  observations  were  available  for 
each  grid. 


Figure  12.  EDAS  analysis  of  500-hPa  heights  and  winds  valid  at  1200  UTC,  27  March  1998.  Solid 
lines  are  geopotential  heights  in  DM;  maximum  wind  vector  =  47.5  m/s. 
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Figure  13.  EDAS  analysis  of  500-hPa  heights  and  winds  valid  at  1200  UTC,  29  March  1998. 
Solid  lines  are  geopotential  heights  in  DM;  maximum  wind  vector  =  54.5  m/s. 


Figure  14.  EDAS  analysis  of  500-hPa  heights  and  winds  valid  at  1200  UTC, 
31  March  1998.  Solid  lines  are  geopotential  heights  in  DM; 
maximum  wind  vector  =  59.0  m/s. 
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The  statistical  parameters  produced  are  the  mean  absolute  error,  root-mean  square  error, 
correlation  coefficient,  standard  deviation,  and  bias;  for  brevity,  only  the  mean  absolute  error  is 
extensively  discussed  in  this  section. 

Table  2  shows  the  BFM  mean  absolute  windspeed  errors  for  all  seven  forecasts,  classified  by 
grid  region  and  height- above-ground  level  with  the  “/”  symbol  separating  the  12-  and  24-h 
forecast  errors.  Mean  observed  surface  windspeeds  in  the  study  period  ranged  from  2.8  m/s 
(Washington  State)  to  7.3  m/s  (Ft.  Hood,  TX).  For  each  of  the  five  grids  that  had  average 
observed  speeds  greater  than  5  m/s,  the  BFM  exhibited  a  slight  underforecast  bias  of  generally  1 
to  2  m/s  (not  shown).  Above  the  surface,  table  2  shows  a  distinct  increase  in  the  mean  absolute 
windspeed  errors  in  proximity  to  the  locations  of  the  lower  and  upper  level  jets,  associated  with 
the  larger  observed  windspeed  values  at  these  levels.  To  examine  these  errors  more  closely,  the 
nearest  NOGAPS  grid  points  to  each  radiosonde  location  were  checked  for  the  New  York  study 
area. 

Tables  3,  4,  and  5  compare  the  NOGAPS  mean  absolute  errors  for  both  windspeed  and 
temperature  at  the  grid  points  closest  to  each  of  the  radiosonde  sites:  Brookhaven,  NY  (OKX), 
Chatham,  NY  (CHH),  and  Albany,  NY  (ALB).  The  975-hPa  level  at  site  OKX  is  examined  due 
to  its  proximity  throughout  the  period  to  a  low-level  windspeed  maxima  and  a  temperature 
inversion;  the  925-hPa  level  at  sites  CHH  and  ALB  was  examined  for  similar  reasons.  During 
much  of  this  period,  a  strong  southwesterly  flow  between  975  hPa  and  925  hPa  was  present  over 
portions  of  the  Northeast  (figure  15).  By  comparing  the  values  in  table  2  with  those  in  tables  3  to 
5,  it  appears  that  the  BFM  500-magl  windspeed  errors  over  the  grid  are  related  to  NOGAPS  low- 
level  windspeed  errors  (used  for  nudging).  However,  a  few  case  studies  of  the  BFM  since  this 
work  have  revealed  that  speed  errors  in  the  vicinity  of  jet  streaks  might  also  be  linked  to  use  of 
too  strong  a  nudging  coefficient  (3  to  5  x  10'"^)  in  the  prognostic  momentum  equations.  Perhaps 
in  the  presence  of  large-scale  horizontal  windspeed  gradients,  the  BFM  may  be  “double 
counting”  the  effect  of  the  horizontal  advection  by  using  a  stronger  nudging  term.  Note  that  this 
idea  still  needs  to  be  investigated  more  thoroughly.  The  nudging  coefficient  value  of  3  x  lO'"^, 
used  in  this  study,  may  explain  some  of  the  larger  BFM  500-hPa  wind  errors  noted  over  several 
of  the  study  areas. 

Tables  6  to  9  compare  mean  BFM  absolute  errors  for  several  parameters,  showing  that  maximum 
errors  tend  to  occur  near  500  magi,  1500  magi,  and  500  hPa  for  most  study  areas.  The  BFM 
forecast  errors  at  lower  levels  may  be  attributed,  to  some  degree,  to  the  relatively  coarse  vertical 
resolution  of  the  NOGAPS,  small  spatial  and  temporal  errors  in  NOGAPS  forecasted  synoptic 
features,  and  nesting  scale-separation  issues.  The  presence  of  strong  low-level  jets  (figure  15)  or 
shallow,  cold  air  polar  surges  offer  other  problematic  situations.  For  example,  within  the  Denver 
grid  study  area,  large  500-magl  wind  errors  were  generated  on  two  occasions  when  observed 
low-level  10  m/s  northeasterly  upslope  flow  was  observed  near  Denver.  On  each  occasion,  both 
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models  (BFM  and  NOGAPS)  forecasted  winds  of  nearly  10  m/s  from  the  southwest. 


At  higher  levels  in  the  atmosphere,  larger  windspeed  errors  are  related  to  errors  in  the  NOGAPS 
forecasts,  in  addition  to  the  strength  of  the  nudging  coefficient  applied.  As  discussed  earlier,  the 
subsequent  use  of  a  weaker  nudging  coefficient  over  Denver  resulted  in  superior  500-hPa  wind 
forecasts  over  those  shown  in  this  report.  Table  10  lists  the  composite  mean  absolute  error  and 
correlation  coefficient  statistics  for  all  regions  by  parameter  and  level,  with  overall  mean 
absolute  errors  lowest  and  correlation  coefficients  highest  at  the  surface  (due  to  the  stronger 
influence  of  terrain  and  surface  feature  irregularities). 


Table  2.  Mean  absolute  BFM  windspeed  errors  (m/s)  by  grid  region  from  27-30  March  1998.  Errors  are  for  12/24-h  forecasts,  and 
SFC  refers  to  10-magl  level. 


Level 

New 

York 

Southern 

California 

Southern 

NM 

Louisiana 

Wash. 

State 

Chicago 

Denver 

Florida 

Ft.  Hood 

TX 

SFC 

1.68/1.72 

2.31/2.54 

2.70/2.67 

1.48/1.69 

1.91/1.93 

1.69/1.96 

3.01/2.60 

1.67/2.25 

2.18/2.51 

500  magi 

3.34/3.61 

3.19/3.19 

3.45/3.19 

2.25/2.72 

2.10/2.07 

3.21/3.49 

3.78/4.62 

2.23/3.29 

4.92/5.32 

1500  magi 

2.97/2.90 

3.21/2.79 

2.72/2.96 

3.16/2.94 

1.49/1.73 

2.73/2.67 

1.59/2.79 

1.67/2.57 

3.51/3.75 

700  hPa 

2.41/1.99 

3.30/3.49 

2.89/3.17 

2.48/2.92 

1.53/2.31 

3.62/2.57 

2.73/3.79 

1.33/1.69 

2.76/2.56 

500  hPa 

2.90/2.63 

2.77/2.78 

3.22/2.85 

2.43/3.80 

2.69/2.93 

2.63/3.40 

4.41/3.88 

2.93/2.87 

5.36/7.22 
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Table  3.  Mean  absolute  errors  of  NOGAPS  point 

41°N,  73  °W  compared  with  OKX  sounding 
in  New  York  grid  area. 


Level 

Wind  speed  (m/s) 

Temperature  (°C) 

975  hPa 

6.35/5.48 

4.57/5.23 

850  hPa 

1.99/0.90 

1.83/2.49 

700  hPa 

2.94/2.50 

1.21/2.08 

500  hPa 

2.09/2.35 

0.87/1.08 

Table  4.  Mean  absolute  errors  of  NOGAPS  point 

42°N,  70°W  compared  with  CHH  sounding 
in  New  York  grid  area. 


Level 

Wind  speed  (m/s) 

Temperature  (°C) 

925  hPa 

2.40/2.73 

2.65/3.85 

850  hPa 

2.61/2.22 

2.23/2.13 

700  hPa 

2.36/2.74 

0.62/1.03 

500  hPa 

3.52/1.75 

0.95/1.17 

Table  5.  Mean  absolute  errors  of  NOGAPS  point 

43  °N,  74° W  compared  with  ALB  sounding 
in  New  York  grid  area. 


Level 

Wind  speed  (m/s) 

Temperature  (°C) 

925  hPa 

3.31/3.76 

4.04/4.44 

850  hPa 

4.25/5.94 

2.16/2.69 

700  hPa 

3.17/2.61 

1.09/1.23 

500  hPa 

3.04/3.15 

1.05/1.43 
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Figure  15.  EDAS  analysis  of  925-hPa  wind  vectors  valid  at  0000  UTC,  29  March  1998.  Maximum 
wind  vector  =  26.0  tn/s  over  southeast  Massachusetts. 
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Table  6.  Mean  absolute  BFM  temperature  errors  (°C)  by  grid  region  from  27-30  March  1998.  Errors  are  for  12/24  h  forecasts, 
and  SFC  refers  to  2-magl  level. 


Level 

New 

York 

Southern 

California 

Southern 

NM 

Louisiana 

Wash 

State 

Chicago 

Denver 

Florida 

Ft.  Hood, 
TX 

SFC 

2.88/2.83 

2.99/2.84 

2.84/3.47 

1.47/1.33 

2.61/2.29 

2.84/3.49 

2.79/2.99 

2.79/2.86 

1.68/1.79 

500  magi 

3.52/5.16 

1.05/1.11 

2.91/3.63 

1.23/1.10 

1.18/1.11 

2.57/2.96 

2.56/2.29 

0.95/0.77 

1.91/2.07 

1 500  magi 

1.59/2.11 

1.16/1.51 

2.05/2.50 

0.86/1.19 

0.82/0.68 

0.70/0.76 

1.76/1.73 

1.36/0.91 

1.71/2.43 

700  hPa 

1.11/1.07 

0.93/1.14 

2.18/2.40 

1.63/1.73 

1.16/1.45 

1.01/1.42 

1.82/1.55 

1.08/1.15 

1.75/2.03 

500  hPa 

0.77/1.11 

1.38/1.87 

1.26/1.69 

0.65/0.89 

1.50/1.50 

0.81/0.84 

1.84/2.69 

0.60/0.58 

0.97/1.02 

Table  7.  Mean  absolute  BFM  u-component  wind  errors  (m/s)  by  grid  region  from  27-30  March  1998.  Errors  are  for  12/24-h 
forecasts,  and  SFC  refers  to  10-magl  level. 


Level 

New 

York 

Southern 

California 

Southern 

NM 

Louisiana 

Wash 

State 

Chicago 

Denver 

Florida 

Ft.  Hood, 
TX 

SFC 

2.18/2.28 

3.07/3.26 

3.08/3.34 

1.13/1.17 

2.77/2.72 

1.85/2.04 

4.17/4.14 

1.84/2.38 

1.80/1.89 

500  magi 

2.60/2.58 

2.92/2.84 

3.43/3.28 

1.41/1.72 

2.75/3.10 

3.07/4.53 

6.16/7.94 

2.80/3.63 

2.47/3.04 

1 500  magi 

2.91/2.70 

2.99/2.80 

2.74/3.30 

1.99/2.21 

1.95/2.87 

2.75/3.72 

3.55/5.14 

1.84/2.92 

3.12/2.80 

700  hPa 

2.28/1.89 

2.64/2.19 

3.09/3.56 

1.73/2.01 

2.15/2.97 

4.09/4.26 

4.56/6.37 

1.04/1.47 

3.25/3.05 

500  hPa 

3.19/2.92 

3.17/3.24 

4.69/4.64 

2.36/2.93 

2.15/3.09 

2.82/2.99 

5.45/5.53 

3.52/3.17 

4.14/4.25 
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Table  8.  Mean  absolute  BFM  v-component  wind  errors  (m/s)  by  grid  region  from  27-30  March  1998.  Errors  are  for  12/24-h 
forecasts,  and  SFC  refers  to  10-magl  level. 


Level 

New 

York 

Southern 

California 

Southern 

NM 

Louisiana 

Wash 

State 

Chicago 

Denver 

Florida 

Ft.  Hood, 
TX 

SFC 

1.75/1.52 

2.57/2.47 

3.21/3.26 

1.55/1.74 

2.12/1.85 

1.45/2.34 

3.57/3.24 

1.22/1.30 

2.52/2.71 

500  magi 

3.44/3.37 

3.33/3.47 

2.77/3.36 

2.23/2.62 

2.68/2.49 

3.17/3.63 

5.43/4.98 

1.36/1.60 

5.39/5.93 

1 500  magi 

2.65/2.86 

3.41/3.77 

3.45/3.15 

3.25/3.06 

1.48/1.61 

2.81/3.33 

1.96/2.01 

1.28/1.44 

4.44/4.71 

700  hPa 

2.91/2.55 

4.04/4.46 

3.74/3.66 

2.52/3.13 

1.09/1.61 

2.59/3.03 

1.89/1.98 

1.34/1.97 

2.85/2.92 

500  hPa 

2.59/2.29 

4.11/4.18 

3.84/2.95 

2.48/4.23 

3.24/2.64 

2.26/3.25 

1.46/2.76 

2.44/2.88 

4.93/8.88 

Table  9.  Mean  absolute  BFM  dewpoint  errors  (°C)  by  grid  region  from  27-30  March  1998.  Errors  are  for  12/24-h  forecasts,  and 
SFC  refers  to  2-magl  level. 


Level 

New 

York 

Southern 

California 

Southern 

NM 

Louisiana 

Wash 

State 

Chicago 

Denver 

Florida 

Ft.  Hood, 
TX 

SFC 

2.52/1.96 

2.29/2.27 

2.41/2.87 

1.56/1.49 

2.73/2.55 

1.93/1.84 

2.22/2.07 

1.08/1.14 

2.43/2.98 

500  magi 

2.12/2.19 

1.78/1.86 

1.99/2.21 

1.45/1.43 

1.45/1.10 

2.41/2.21 

2.43/2.61 

1.60/1.28 

2.36/2.06 

1 500  magi 

1.74/2.13 

2.23/2.60 

2.42/2.47 

3.11/2.80 

1.87/1.90 

2.24/2.44 

1.99/1.96 

4.11/3.48 

3.60/4.03 

700  hPa 

2.58/2.31 

3.60/3.35 

2.43/2.31 

1.80/2.13 

2.31/3.58 

2.33/2.94 

2.27/2.16 

1.34/1.40 

2.92/3.62 

500  hPa 

1.98/1.41 

2.18/3.17 

3.19/3.16 

0.69/1.42 

2.78/3.06 

2.05/2.33 

2.22/2.33 

0.60/0.58 

0.97/1.02 

28 


Table  10.  Composite  mean  absolute  error/correlation  coefficient  over  all  regions  and 
times  for  selected  parameters  and  levels. 


Level 

Wind  Spd 
(m/s) 

U-comp 

(m/s) 

V-comp 

(m/s) 

Temp 

rc) 

Dew  pt 
(°C) 

SFC 

1.82/0.55 

1.98/0.55 

1.96/0.63 

2.35/0.83 

1.78/0.72 

250  magi 

2.28/0.64 

2.44/0.70 

2.34/0.80 

1.60/0.90 

1.41/0.78 

500  magi 

2.32/0.71 

2.49/0.72 

2.55/0.84 

1.47/0.90 

1.38/0.75 

700  hPa 

2.10/0.81 

2.25/0.83 

2.42/0.91 

1.16/0.89 

1.84/0.67 
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6.  Possible  Improvements  and  Considerations 


Several  features  of  the  eurrent  BFM  eould  be  examined  for  improvements  and  modifieations — 
sueh  as  the  use  of  the  Mesoseale  Modeling  System  V  (MM5)  of  the  Pennsylvania  State 
University/National  Center  for  Atmospherie  Researeh  (Grell  et  ah,  1994)  for  obtaining  nudging 
input  files.  AFWA  is  now  operationally  running  the  MM5  at  a  45-km  horizontal  resolution, 
offering  the  potential  to  reduee  the  seale-separation  nesting  ratio  for  nudging  to 
4.5-to-l,  a  signifieant  improvement  from  the  eurrent  eonfiguration  ratio  using  NOGAPS. 

Other  possible  improvements  might  involve  a  better  treatment  of  elouds  and  radiation,  an 
elevation-dependent  weighting  for  determining  the  surfaee  nudging  eoeffieients,  a  higher 
resolution  and  heterogeneous  surfaee  state  input,  and  a  eumulus  parameterization  option  (Sun 
and  Haines,  1996). 


7.  Summary  and  Discussion 


The  BFM  is  a  useful  short-range,  operational  foreeasting  and  planning  tool  for  the  taetieal  Army, 
despite  its  limitations  (i.e.,  loealized  summer  eonveetion,  nonhydrostatie  foreing,  sharp  frontal 
boundaries,  ete.).  The  use  of  Newtonian  relaxation  provides  an  ability  to  fully  assimilate  foreeast 
tendeneies  from  a  larger-seale  numerieal  model.  In  benign  synoptie  regimes,  sueeessful 
simulations  of  loeal  boundary  layer  flows  foreed  by  surfaee  inhomogeneities  in  terrain  and 
radiation  are  possible.  If,  for  example,  the  AFWA  MM5  were  to  be  used  to  obtain  BFM  nudging 
guidanee  (replaeing  NOGAPS),  the  horizontal  resolution  nesting  ratio  eould  be  redueed  to  about 
4-to-l  or  greater.  This  would  reduee  the  seale  separation  and  resolution  eonsisteney  problems 
diseussed  by  Warner  et  al.  (1997).  Nonetheless,  previous  results  have  shown  that  the  present 
BFM  ean  still  provide  useful  foreeasts  in  the  lower  boundary  layer  (Haines,  1999),  whieh  ean 
guide  taetieal  weather  deeision  aids  and  other  weather-dependent  eomponents  of  the  Intelligent 
Planning  of  the  Battlefield  proeess. 
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Acronyms  and  Abbreviations 


ADI 

alternating  direction  implicit 

AFWA 

Air  Eorce  Weather  Agency 

AI 

artificial  intelligence 

ALB 

Albany  radiosonde 

ASP 

Atmospheric  Sounding  Program 

BFM 

Battlescale  Eorecast  Model 

CAAM 

Computer  Assisted  Artillery  Meteorology 

CHH 

Chatham  radiosonde 

CFL 

Courant-Eriedrich-Eewy 

EDAS 

ETA  Data  Analysis  System 

FCMM 

forecast  computer  meteorological  messages 

GMDB 

gridded  meteorological  database 

GTRAJ3 

General  Trajectory  Model  3 

GUI 

graphical  user  interface 

IMETS 

Integrated  Meteorological  System 

IWEDA 

Integrated  Weather  Effects  Decision  Aid 

EUT 

Eimited  User  Test 

Met 

meteorological 

MM5 

Mesoscale  Modeling  System  V 

MRMD 

mean  radial  miss  distance 

NOGAPS 

Naval  Operational  Global  Atmospheric  Prediction  System 

NWP 

numerical  weather  prediction 
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OKX 

Brookhaven  radiosonde 

OWS 

Operational  Weather  Squadron 

PBL 

planetary  boundary  layer 

RCMM 

radiosonde  computer  meteorological  messages 

RDAP 

Reliability  Determination  and  Assurance  Program 

SAD ARM 

Sense  and  Destroy  Armor 

TDAs 

Tactical  Decision  Aids 

T-VSAT 

Tactical- Very  Small  Aperture  Terminal 

UTC 

universal  time  coordinated 

XMR 

Cape  Canaveral  radiosonde 

YPG 

Yuma  Proving  Ground 

4DDA 

four-dimensional  data  assimilation 
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